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PREFACE 


This report is the second of the two volume documentation for 
the bearing analysis computer program SASHBEAN. This volume, 
Volume II, provides the details of the underlying mathematical 
formulation, analysis, and solution algorithms used for this 
computer program. A separate volume. Volume I, provides the 
detailed instructions required to successfully install and 
effectively use the software for the design and analysis of 
single row, angular contact, spherical roller bearings. 

All efforts involved in the development of this software and its 
documentation were performed by McGill Manufacturing, Emerson 
Power Transmission Corporation. This work was done as part of 
the Advanced Rotorcraft Transmission (ART) Program to advance 
the state-of-the-art in helicopter transmissions. The ART 
program was funded by the U.S. Army Aviation Systems Command 
(AVSCOM) and managed cooperatively by the AVSCOM Propulsion 
Directorate and the NASA Mechanical Systems Technology Branch, 
both located at the NASA Lewis Research Center, Cleveland, Ohio. 
This work was done under a sub-contract to Sikorsky Aircraft 
Division of United Technologies Corporation, the prime 
contractor, under NASA contract NAS3-25423. 

Technical direction for this project was provided by Sikorsky 
Aircraft's representatives Mr. C.H. Keller, Jr. and 
Mr. J.G. Kish, the Task Manager of the project. The 
government's technical representatives for this work were 
Dr. R.C. Bill, ART Program Manager and Mr. T.L. Krantz, Project 
Manager for the Sikorsky ART contract. 

The activities performed at McGill Manufacturing were directed 
by Mr. D.M. Michaels, Project Manager for the sub-contracted 
project. Analytical and technical support was provided by 
Mr. J.S. Porter, Mr. R.H. Barber, Mr. C.A. Kruse, 

Mr. G.A. Satkamp, Mr. A.K. Aggarwal and Mr. W.D. Nutt. Drawing 
and drafting aid were provided by Mr. D. Wisch and 
Mr. T. Peterson. Typing and word processing were done by 
Ms. C. Dodrill and Ms. B. Richards. 
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LIST OP SYMBOLS 


Alphabets 


a 

A 

b 

c 

C 

CF 

d 

dc 

D 

E 

E 

e 

F 

f 

GM 

h 

H 



i 


j 

1 

K 


Semi -major axis of the contact ellipse (in) 
Dimensionless material parameter for EHD analysis 
Semi-minor axis of the contact ellipse (in) 

Specific heat of a material (Btu/lbm. ®F) 

Initial clearance between a roller -raceway laminae (in) 
Centrifugal force acting on each roller (lb) 

Distance (in) 

Diametral clearance (in) 

Roller diameter, maximum (in) 

Young's modulus of elasticity (psi) 

Equivalent modulus of elasticity (psi) 

Poisson's ratio of a material 

External load on the bearing (radial /axial) (lb) 
Friction force at a sliding contact (lb) 

Gyroscopic moment acting on each roller (in- lb) 

E HP film thickness, minimum (in) 

Rate of heat generation (Btu/hr) 

Moment of Inertia of a roller about it's 
longitudinal axis (in- lb sec ) 

Moment of Inertia of a roller about it's transverse 
axis (in- lb sec^) 

Ring ID; i = 1 for inner; i = 2 for outer 
Roller ID; 1 s j s Z 

Lamina ID; 1 s 1 s F 

Material constant in the load- deformation relationship 
for a lamina contact 
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L 


L-10 = 


LP 

M 

N 

n = 


P 

P 


P 


e 



Q 

q 


r 


R 

S 

t 

T 

U 

u 

V 

W 

w 

y 

z 


Roller length (in) 

L-10 life of a raceway (million revolutions) 

Lvibricant parameter of the liabricant (sec) 

Total pitching moment acting on a roller (in- lb) 
Rotational speed of the bearing (rpm) 

Total number of loaded laminae at a raceway 
Mass density of a material (Ibm/in^) 

Resultant roller load (radial /axial) (lb) 

"Equivalent" roller load at a roller- raceway 
contact (lb) 

"Mean" raceway load for all rollers (lb) 

Basic dynamic capacity of a raceway (lb) 

Normal load at a lamina contact (lb) 

Normal load per unit length for a line contact (lb) 
Radius of curvature of the contacting surface (positive 
for a convex surface, negative for a concave surface) 
Crown radius of roller or raceway 
Contact stress (psi) 

Time (sec) 

Viscous torque of the lubricating fluid (in- lb) 

Mean entrainment velocity (in/sec) 

Dimensionless speed parameter for EHD analysis 
Tangential velocity of a rotating point (in/sec) 

Total width of cage rails (in) 

Width of each lamina (in) 

Dimensionless load parameter for EHD analysis 
Total number of rollers in the bearing 


5 



Greek Letters 


<t> 

T 

T 

5 
$ 

Of 

6 

! 

5 

e 

6 

"k 


ff 

/ 

Q 

M 


Roller azinvuth angle (radians) , 0 s 0 s 2ir 
Shear stress (psi) 

Total n\jinber of lamina for a roller 
Normal approach at a line contact (in) 

Ring deflection, relative (in) 

Roller angular displacement (pitching) (in/ in) 
Initial contact angle of the bearing (radians) 
Operating contact angle of a lamina at 
the raceway (radians) 

D.CosS/Pt (bearing geometry parameter) 

Roller deflection (radial /axial) (in) 

Kinematic viscosity of the lubricant (cStokes) 

(@ operating temperature and atmospheric pressure) 
Absolute viscosity of the lubricant (Ib.sec/in ) 

(@ operating temperature and atmospheric pressure) 
Surface finish, RMS (in) 

Liibrication film parameter from EHD analysis 
Angular velocity of rotation (rad/sec) 

Coefficient of sliding friction 


Special Characters 
ae = Weibull slope 

JE = Material factor in the load- deformation 

relationship for a line contact 
c = Pressure -Viscosity coeff. of the lubricant (in /lb) 
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Pitch diameter of the bearing (in) 

Basic dynamic capacity reduction factor 
in fatigue life equation 

Empirical factor based on bearing type and it's 
l\ibrication system. 

Indicates functional relationship 


Subscripts 


a = 

brg = 
c 

eg 

cr 

eff = 
fl 

i = 

j 

1 

m = 
max = 
mean = 
min = 
o 
r 

tr 


Axial direction 

For the coit^lete bearing 

Composite 

For cage /retainer 

For cage rails 

Effective 

For lubricating fluid 

Ring niunber; i = 1 for inner, i = 2 for outer 

Roller nximber; 1 s j s Z 

Lamina ntimber; 1 s 1 s r 

Roller mid-plane 

Maximum value 

Mean value 

Minimum value 

Any roller 

Radial direction 

Transition point 
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1.0 INTRODUCTION 


Spherical roller bearings, known for their high load carrying 
capacity along with their ability to perform under conditions of 
misalignments and shaft deflections, have conventionally been 
restricted to "low" to "moderate” speed applications. Speeds to 
the order of 5000 rpm or 0.25 million DN have typically [17] been 
considered as the upper limits for this class of bearings. 

The higher coefficient of friction exhibited by a spherical roller 
bearing, as compared to that of an "equivalent” size cylindrical 
roller bearing, has been one of the main deterrents for the 
former's use in high speed applications. In these bearings, the 
close conformity of roller and raceway spherical crowns, while 
lending the bearing its self-aligning and high load capacity, also 
results in relative sliding at he concentrated contacts, resulting 
in the higher overall coefficients off rolling friction for the 
bearing. 

New concepts in materials, manufacturing techniques, lubricants, 
component design and design tools are causing revolutions in 
bearing performance, and spherical roller bearings are no 
exception. Spherical roller bearings arc now being designed and 
developed for high load and/or high speed applications including 
aerospace applications. 

One such design has been developed by McGill Manufacturing, under 
a sub-contract from Sikorsky Aircraft Company, as part of an 
ARMY/NASA sponsored ART project. The bearing, approaching 1.15 
million DN, has been tested under full load and full speed 
conditions . 

The trend towards lightweight, high speed, and high performance 
applications, with increased requirements of reliability and 
safety, has also placed a great deal of emphasis on the ability to 
accurately analyze and predict the performance of any suggested 
design. Such rigorous bearing analysis is no longer restricted by 
the availability of main frame and mini computers. With the 
phenomenal development in the power and speed of desk top personal 
computers, such detailed analysis is now practical on personal 
computers . 

The McGill computer program, SASHBEAN, based on the mathematical 
formulation described in this report, provides a sophisticated 
analytical tool to design, analyze, and predict the operating 
characteristics of single row, angular contact (including zero 
degree contact angle) spherical roller bearings under high speed 
conditions. 
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2.0 FORMULATION OF MATHEMATICAL MODEL 


To simulate the dynamic performance characteristics of a bearing, 
the mechanics of internal motions, load and stress distributions, 
lubricant flow, sliding friction, and heat generation have to be 
modeled mathematically and solved for the given parameters of 
bearing geometry, material, empirical factors, and external load 
and speed environment. 


When the bearings involved are of special design, incorporating 
non-convent ional materials, and operating under ertreme ^ 
conditions, these mathematical models, with few simplifying 
assumptions, no longer allow hand calculations. Computer programs 
running on high speed digital computers and employing efficient 
numerical solution techniques are often required to solve such 
complex mathematical models. 


The mathematical 
program SASHBEAN 


model that forms the basis of the computer 
is described in detail in the following sections. 
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3.0 LOAD DISTRIBUTION ANALYSIS 


As in any bearing analysis, a major effort involves modeling and 
so3_ving for the bearing^ s load distribution among its rolling 
elements for the given load and speed environment. The loading of 
the rings at many rolling elements poses a statically 
indeterminate problem and is often complicated to solve. The 
problem is further complicated when the high speed dynamic loads 
on the rolling elements, namely centrifugal forces and gyroscopic 
moments, are not negligible and are fully considered. 

A method, referred to here as the "LAMINA" method, has been 
employed to solve for the load distribution in the bearing ^ under 
the combined environment of externally applied loads and high 
speed dynamic forces. In this technique the roller is divided 
into a number of slices (laminae), the slicing planes being normal 
to the roller axis. Similarly each raceway is also considered to 
be made up of an equal number of laminae. As the deformation at 
any concentrated contact is very small, it is further assumed that 
the inter-facial shear between a loaded lamina and an adjoining 
non-loaded lamina is negligible and that the loaded lamina deforms 
independently of the non-loaded one under normal loading. 

3.1 COORDINATE SYSTEM AND SIGN CONVENTION 

Consider a single row, (non-zero) angular contact spherical roller 
bearing as shown in Figure 1. Let he outer ring be stationary and 
the inner ring rotating for this formulation. In a real life 
situation, if the outer was rotating or both inner and outer were 
rotating, the same formulation will hold as the roller dynamic 
loads are estimated based on the actual rotational speed of the 
two rings . 

The rollers are numbered 1 through Z sequentially in the clockwise 
direction and equally spaced, as shown. The azimuth (angular) 
location of roller #1 is considered as 0.0 degrees and the same 
increases in the clockwise direction with the roller number. 
Therefore, the azimuth angle, (j>, of any roller, j, is given by, 

= 2 T(j-l)/z, 1 < j < Z (3.1) 

Let the radial load on the bearing be acting along the Y-axis, 
with roller #1 being directly under the load as shown in Figure 1. 
The relative radial deflection of the rings is thus positive along 
the positive direction of the radial load. 

An axial load, applied either at he inner or outer ring, is 
positive when the two rings are pressed into each other. Accord- 
ingly, the relative axial deflection of rings is positive when 
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this deflection causes the rings to move ^ially towards each 
other. An axial deflection causing the rxngs to move away from 
each other, from the zero end-play position, is thus in the 
negative direction. A positive axial load, acting along the 
X-axis is also shown in Figure 1. 

Due to the self -aligning ability of a spherical roller bearing, 
any external moment loading, if present, is not considered. 
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3 .2 CONTACT LOAD - DEFORMATION RELATIONSHIP 


Por a cylindrical body of finite length, L, when pressed on to a 
plane surfaced body of infinite extension, the normal approach, 

5, between the axis of the cylinder and a distant point in the 
supporting body is approximately given by, 

d = 4.36E-7 (3.2) 


where JE is a factor based on the materials of the two contacting 
bodies q is the normal force per unit length of the contact- 

ing cylinder. Appendix E describes in detail the computations 
for the material factors for different contacting materials. 


If the effective length, Lg^^, of the roller is s\ibdivided 
(sliced) in r laminae (slices) , each of width w, then equation 
(3.2) can be rewritten as. 


or 


a = 4.36E-7 [JE2’*^q*^ (rw) '^] 

Q = 


{gl.llw.89) 


(Kr*^^) 


(3.3) 


where K = (4 . 36E- 7JE^ * *^) ^ * ^^ is another constant and Q = q.w is 
the total normal force at a lamina contact. 


3.3 initial CLEARANCES BETWEEN ROLLER - RACEWAY LAMINAE 


Due to the difference in the radii of curvatures of roller and 
raceway crown profiles, there exists varying initial clearances 
between the roller and raceway laminae. When the bearing is 
loaded, the contact of a roller lamina, 1, with that of a raceway 
lamina takes place only after this initial clearance is removed. 
Figure 2 shows the relative position of a roller, j, with respect 
to the two raceways when the bearing is held together with zero 
radial and axial deflection of the rings. As the roller is sub- 
divided into r lamina, each of width w, the initial clearances at 
any lamina, 1, can be determined from geometry as follows. 


Ci,i = Rq- (R o^-d^^)^- [R^- (Ri^-dj^^)^] , 1 s 1 s r & 

Where di = w(l-lj^) is the distance of the lamina, 1, 
roller mid-plane? Ij^ being the roller's mid-lamina. 


i=l,2 (3.4) 

from the 


Having established the contact load- deformation relationship and 
the initial clearances present between the roller and raceway 
laminae, the load distribution analysis then reduces to determin- 
ing the relative axial and radial deflections of the rings when 
in equilibrixam under the externally applied loads and the inter- 
nal forces. For these ring deflections, each roller has to be in 
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operating equilibrium under the contact forces at raceways and 
the dynamic loads . 


AS the loading of the rings at the many rolling element points 

^ses a statically indeterminate probl^, mCSl^'^arTteratiJe 
the inclusion of the dynaitixc loads in the . . 

so!u“on sJh^e has bsS deployed to solve this load dxstributton 

problem. 



FIGURE 2 


13 




3.4 THE ITERATIVE SOLUTION 


Ttis bGaring irdLn^s slxq 3 .ssic['n©d. s. sm3.Xl but known sst of nslstivs 
radial and. axial dsflsctions. Let these deflections be and 
respectively. For the given ring deflections, and we then 
determine the radial, axial, and angular deflectxons of each 
roller by satisfying its equilibrium under radial, axial forces 
and pitching (misaligning) moments. This is also done using an 
iterative scheme as described below. 

A roller, j, at azimuth 0 is assigned a given set of axial, 
radial and angular displacements as shown in Figure 3 . Let a, 

€ , e_ be this roller's angular, axial, and radial displace- 

mints respectively. From geometry, for the given ring and roller 
displacements, we can now write expressions for the normal ap- 
proach at each lamina as represented below, 

5. . , = / (Geometry, 4-3, $j.Cos0, €-., €^, a) - 1 (3.5) 

where C- t is the initial clearance at this lamina contact. 

X » X 
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A positive and non-zero 6 . . i then indicates a loaded lamina 
contact If n^ • and n, the number of loaded laminae, for 

the roller, the i^ir and outer respectively, we can set up 

the roller equilibrium equations in the axial, radial, and angu- 
lar directions as follows: 

For the equilibrium of this roller in the axial direction, 

? (3.6) 

I ^ J • 


w n 


•a, 1,3 


Kfid ^ 1=1 

■^ / J 


For the equilibrium of this roller in the radial direction. 


w ' n 


•r,i ,3 


Kn 


i,: 


S 5.i 
1=1 


1.11 


SinSi 


- J-i n 1 -i 1 

.11 10»1 -L/J/J- 


(3.7) 


For the equilibrium of this roller in the angular direction, 

,89 


w' 


M 


1/1 


n 
r 5q 


Kn 


i/1 


.11 


1 . 11 , 


(3.8) 


Where S- • ^ is the contact angle at the raceway, i, roller, j, 
and laming 1 . 

With the contact forces and moments now known for roller, j, its 
(dynamic) equilibrium condition is then established as follows. 

a. Check equilibrium of this roller in axial direction: 

l^^a2,j ■ ^a,l,j^l > Allowed Tolerance (3.9) 

THEN adjust this roller's axial deflection and start over for 
this roller. 

b. ELSE Check e<guilibrium of this roller in radial direction: 

jp I (p - p , .) - CFj > Allowed Tolerance (3.10) 

THEN adjust this roller's radial deflection and start over for 
this roller. 


c. ELSE check equilibriiom of this roller's pitching moments: 

IF 1 (Mo • + M, ^) - GMl > Allowed Tolerance (3.11) 

THEN adjust this roller's angular displacement and start over for 
this roller. 

ELSE the roller, j, is found to be in equilibrixam for the given 
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axial, radial and angular roller displacements. 

This process is then repeated for all rollers in the bearing and 
their equilibrium positions determined. In actual programming, 
taking advantage of the bearing symmetry about a plane through 
the bearing axis, only a part of the actual number of rollers are 
solved . 

With the equilibrixam forces for all the rollers in the bearing 
now known for the assigned ring deflections, the overall bearing 
equilibrium equations against the externally applied radial and 
axial loads are then set up as follows, 

(a) Check bearing equilibrium in the radial direction: 


IF 


E P 

j=l 


r,l, j 


Cos0j 



> Allowed Tolerance 


(3.12) 


THEN adjust rings ' relative radial displacement and start over a 
coir 5 >lete new iteration. 

(b) ELSE check bearing equilibrixim in the axial direction: 


IF 


Z 

E P 

j=l 


a, 2, 


j. 


> Allowed Tolerance 


(3.13) . 


THEN adjust rings ' relative axial displacement and start over a 
complete new iteration. 


ELSE equilibrium of the bearing is established under the given 
loading environment for these ring deflections . 


The normal loads and deflections at each lamina contact are thus 
known for this (equilibrivim) state of the bearing. 


A graphical representation (flow chart) of the iterative scheme 
for this load distribution analysis is shown in Figure 4. 
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4.0 CONTACT STRESS ANALYSIS 


4.0 MAXIMUM AND MEAN CONTACT STRESSES 


Having determined the normal forces at each lamina contact , 
at inner and outer raceways, the level of contact stresses and 
areas are then estimated. The simplified formulation of Hertz 
contact stress analysis for line contacts has been used. The 
roller and the raceway laminae (slices) are considered as two 
cylinders of equal length (Icimina width), pressed against each 
other under a known normal load. As the contact deformation at 
any lamina contact is very small compared to the overall body 
dimensions, the laminae interface shear is assumed to be neglxgi- 
bly small and thus neglected. A typical arrangement of two 
parallel cylinders of equal lengths pressed against each other 
xinder a normal force is shown in Figure 5. The half width of the 
lamina line contact, b, as shown in Figure 5, is given by. 


^ r 1 




7rw(rQ"-^+rj^‘-^) 


4.1) 


where e and are elastic constants of the two bodies based on 
the respective values of modulus of elasticity Poisson's 

ratio of their materials. r_ and r^ are the radii (with proper 
signs as per sign convention) of tne roller and raceway laminae 
respectively as shown in Figure 5. The maximum and mean contact 
stresses are then given by. 


Smax = (2Q)/(7Twb) (4.2, 

Smean = (Q) / (2wb) = (^*3 

In the computer program, the contact width, maximum and mean 
contact stresses are computed at each lamina contact.. The maxi- 
mum of maximxam and the maximum of mean contact stresses of all 
the loaded lamina of the roller, j, at the inner raceway are 
taken as this roller's maximum and mean contact stresses at the 
inner contact . 

Similarly, the maximiim and mean contact stresses of this roller 
the outer contact are given by the maximum of maximum and the 
maximum of mean lamina stresses at the outer contacts . 


18 





4.2 STRESS rONCENTRATION DUE TO EDGE LOADINQ 


If for a roller, any of its "edge" laminae are in contact with 
the raceway and thus loaded, its (roller's) contact ellipse may 
not be fully contained on it's effective length and is truncated 
at the loaded edge. Due to stress concentration, the level of 
contact pressure at such edge laminae is thus expected to be 
higher than that calculated by the above Hertz analysis. 

A user specified stress concentration factor is thus applied 
(multiplied) to these edge laminae normal loads. Accordingly, 
the calculated Hertz contact stresses at these laminae contacts 
also get adjusted by this stress concentration factor. Factors 
of the order of 1.5 - 2.0 are typically used. 

After applying the stress concentration factors, the maximum of 
Tna Y-iTmim and the maximxim of mean laminae contact stresses for each 
roller both at the inner and outer raceways, are also printed out 
■in the program output file. 


4.3 CONTACT ELLIPSE DIMENSIONS 

Cumulative width of the loaded laminae of a roller at the inner 
contacts gives us the major axis of the roller's contact ellipse 
at this raceway. The width of the line contact at the most 
heavily loaded lamina of this roller at the same raceway gives us 
the minor axis of the contact ellipse. 

The axial distance of the most heavily loaded lamina mid- plane 
from the roller mid-plane, being called as the eccentricity of 
the contact ellipse, is also computed and written out for each 
roller. 

The contact ellipse major axis, minor axis, and eccentricity for 
each roller at the outer raceway contacts are also determined in 
a similar fashion and printed out. 

4.4 MAXIMUM SUB - SURFACE SHEAR STRESS AND DEPTH 


As the material below a concentrated contact is also in a state 
of stress and the rolling contact fatigue failures have been 
known to originate from these subsurface points, the itagnitude 
and depth of this subsurface shear stress is also of importance 
to a bearing analyst . 

By considering the stresses caused by the normal contact load and 
further application of the principles of elasticity theory, Jones 
[11] has presented the expressions for the three principal 
stresses occurring at a point along the Z-axis, any depth below 
the contact surface as shown in Figure 6. 
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Since the surface contact pressure is maxirtn^ along the Z- axis, 
the three principal stresses must also attain their maxim^ 
values along the same axis. The maximum shear stress , 

aiven by half of the maximum difference between any two principal 
stressed The depth of this point (of maximum shear stress) can 
also be determined. 

For simplicity, the graphical method suggested by Jones [11] has 
been used for these estimations. The graph presented in 
ence [11], as shown in Figure 7, gives the maxintum shear st re 
and its depth of occurrence as a function of contact ellipse 
dimensions ratio Jb/a. 

with the contact ellipse dimensions and the maximum contact 
pressure known for each roller- raceway contact, the ^imum s^- 
surface shear stress and its depth, are directly read 
graph of Figure 7, in the computer progr^. For a complete 
discussion on this topic and the derivation of underlying equa- 
tions, the reader is referred to Jones [11] and Harris [9] . 



FIGURE 6 
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MAXIMUM SUB-SURFACE SHEAR STRESS/DEPTH 

FROM REFERENCE [11] 



FIGURE 7 


4.5 ROLLER NORMAL LOADS AND OPERATING CONTACT ANGLE S 

The contact angle of the most heavily loaded lamina of a roller 
at the inner is taken as the operating contact angle of this 
roller at the inner raceway. Similarly, the contact angle of the 
most heavily loaded lamina of this roller at the outer is taken 
as the operating contact angle of this roller at the outer race- 
way. 

The siimmation of the components of a roller lamina loads at the 
inner contacts, along the contact angle of the most heavily 
loaded lamina at this raceway, gives the roller- inner contact 
normal load for this roller. Similarly, the roller-outer normal 
load for this roller is given by the summation of l^ina load 
components along the contact angle of the most heavily loaded 
lamina at the outer. 
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5 . 0 EED ANALYSIS 


Grubin's equation [23] has been used forEHD film thictoess 
estimation at each lamina contact. Thxs xs due to the fact that 
the ASME recommended life adjustment curve for EHD lubrxcatxon xs 
also based on computations using the same equatxon as gxven below. 


where , 


h^in = 1.95R_(U) 


eq 


8 / 11 (^) 8 / 11 - 1/11 


(y) 


(5.1) 


= Minimum thickness of the EHD film (in) 

V = fR.R )/(R*±R N is the "equivalent" radius of the two 
cS3tacting°surflceil The plus sign is for external contacts 
(both surfaces convex) and minus for the xntemal contac 
(surface with larger radius of curvature is concave) 

U = (n u) / (ER ) is the dimensionless speed parameter 

A = (<?E) is the dimensionless material parameter 

y = Q/ (ERgqW) is the dimensionless load parameter 

U = ^(Uq+U^) is the mean entrainment velocity (in/sec) 


E 


fl-e, 

= ^ 

L E^ 


2 . - 2 


^ is the "equivalent" modulus of 
E- -J elasticity (psi) 


Rewriting (5.1) with definitions of various parameters, we get 




an 


= 1 . 95Rgq (“a^^^/^eq^ (Q/ERgqW) 


eq 


(5.2) 


Furthermore, in equation (5.2) the two lubricant properties, 
and g can be combined into one parameter, known as Lubrxcant 
Parameter (LP) , as defined below. 


LP = 10 


lln c 
a V 


(5.3) 


This Parameter combines both pressure and temperature -viscosity 
characteristics of the lubricant and thus contaxns ® 

entire contribution to the formatxon of EHD fxlm. Mobxl, for 
example, publishes this data for their premium lubrxcatxng oxls. 

Lubricant properties in the above formulation are expected at the 
o? contacting aurfaces at the ^ 

As the oil film is very thxn, xt xs belxeved that xt quxckly 
attains the inlet surface temperature. If this surface tempera- 
ture is unknown and cannot be easily determxned, the ox ou 
temperature or the average of the inlet and outlet temperatures 
may be used as a good approximation. 
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6.0 FATIGUE LIFE ESTIMATION 


A method, referred to here as "Equivalent Roller Loads" method, 
has been used for the estimation of raceways' L-10 fatigue lives. 
In this method, the nonuniform eind skewed load distribution along 
the roller length, as e3<pected in high speed spherical roller 
bearings, is accounted for by first calculating what is being 
termed as "Equivalent" roller loads. These are computed from the 
roller's individual lamina loads at a raceway using the product 
law of probability. 


According to this law, the probability of survival of the entire 
raceway is the product of the probabilities of survival of each 
individual lamina (slice) . Using the roller "Equivalent" loads 
at a raceway then allows the use of formulation developed for 
estimating the L-10 fatigue life of a roller- raceway under radial 
load with line contact. Harris [9] describes in detail the 
(jerivations for the equation used and presented here for calcu- 
lating the roller "Equivalent" loads, from the Icnown normal 
lamina loads, as follows: 


e,i,D 


= r'^/9 




2/9 


( 6 . 1 ) 


P being the "equivalent" roller- raceway normal contact load at 
raceway, i, and roller, j,. The fatigue lives of inner and outer 
raceway are then given by, 

L-10 = (P^/Pjj^)^ million revolutions 


Where P-, is the Basic Dynamic capacity of the raceway and Pjjj is 
the "mean" roller load at this raceway. P^, for inner and outer 
raceways are given by. 


( 1 - 0 ) 


29/27 


= 49500(£) 


( 1 + 0 ) 


1/4 


LcosS-l 


2/9 


,D)29/27(l^^^) 7/9(2) -1/4 (6.2a) 


(1*6)29/27 9 2/9 


, ■ 49500(£) ■*-[ (D)29/27(i,^^^) 7/9,2)-1/4 (s.2b) 

' (1-6)^/^ LcosfiJ 


Where equation (6.2a) is for the inner and ecpaation (6.2b) is for 
the outer raceway. £ is a Bearing Dynsimic Capacity reduction 
factor based on the bearing type and is a user ^ input to the 
program. Typically, for angular contact spherical roller bear- 
ings, £ lies between 0.60 and 0.85. 

The "mean" roller- raceway load, P^^^, is the quartic mean of the 
roller "equivalent" loads at this raceway and is given by. 
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— — * 

p = r_ 2 P ■ •‘^1 for the rotating raceway (6.3) 

Lz j=l -• 

1 Z n 

p _ r_ r P . . for the stationary raceway (6.4) 

“‘'i ■ Lz j=i -I 

The L-10 fatigue lives of the raceways are then given by, 

L-lOi = (Pc,i/Vi>^ -• ^ 


The L-10 fatigue life of the complete bearing may then be deter 
mined using the law of probability as follows, 

-ae )-l/ae 


L-10^33,g = (L-I0i‘* + L-IO2 


( 6 . 6 ) 


In the computer program, the overall bearxng life is computed and 
printed out, both before and after applying the lubrication and 
"material life adjustment factors to the computed raceway lives. 

6 , 1 life adjustment factor for lubrication 

In the computer program, the h—j^j^ is estimated at each concen- 
trated lamina contact using the user specified data on lubricant 
properties at the expected surface temperatures. The minim^ of 
the minimum EHD film thicknesses at each raceway is then used for 
calculating the film parameter, /, at that 

rms surface finish of the two surfaces is calculated from their 
individual rms surface finishes as follows: 

i = 1, 2 (6.7) 

The film parameter, /, for a raceway is then given by. 


‘'c.i - 


A “ ^in,i/^c,i 


( 6 . 8 ) 


The ASME recommended curve for life adjustment factors has been 
extended for lower values of / (down to v = 0.1) and used for 
reading the raceway life adjustment factors for l^rication. The 
extension of the original ASME curve for lower values of V is 
Ssid on the collected test and field data from Sikorsky. This 
extended ASME curve has been coded into the computer program 
using ctibic splines interpolation and is shown in Figure 8. 
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LUBRICATION LIFE ADJUSTMENT FACTORS 

(ASME CUFWE EXTENDED FOR LOWER LAMBDAS) 



FIGURE 8 


For 1 s / s 10 the life adjustment factor, L/L-10, is returned 
as per the curve shown in figure 8. For / < 0.1, a value of 0.27 
is returned and for > 10, L/L-10 factor of 3.02 is returned. 

The computer program prints out both the estimated and adjusted 
raceways and bearing lives in the program output file. 


6.2 LIFE ADJUSTMENT FACTOR FOR CONSTRUCTION MATERIAL 

An appropriate life adjustment factor, based on the construction 
material of the bearing components, is supplied by the program 
user. This factor is directly multiplied to the calculated L-10 
lives of the bearing raceways along with the life adjustment 
factors due to lubrication. The adjusted raceway lives then 

used in equation (6.6) to estimate the adjusted bearing life. 
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7.0 INTERNAL MOTIONS AND SPEEDS 

With the operating equilibrium position of each roller in the 
bearing now known, the average cage and roller rotation speeds 
are estimated as described below. 


7.1 CAGE ROTATION SPEED 

Assuming no skidding at the concentrated roller- raceway contacts, 
the cage tangential velocity is taken as the average of the 
tLgential vllocities of the most heavily loaded points ^ the 
"controlling" rollers at the inner and outer contacts. The 
rollers loaded both at the inner and outer raceways 
ered as the rollers "controlling" the cage rotation. The rollers 
loaded only at one raceway are, on the other hand, considered a 
"controlled" (orbited) by the rotating cage. 

Consider the operating position of the roller, j, at the inner 
raceway as shown in Figure 9. Let A be the contact center (most 
heavily loaded point) and r^ be the radial distance of point A 
f?om the bearing axis as shJwn. If is the contact angle of 
roller, j, at point A, then from geometry 


= (AP)CosSi 


(7.1) 


Where AP = OP - OA is also determined from geometry. The tangen- 
^II velocity of the point A, when lying on the inner raceway, is 

then given by, _ o ^ (7.2) 


vi = 


Similarly for the most heavily loaded point A at the outer r^ce- 
way contact for this roller, as shown in Figure 10, we can write. 


and 


T2 = (AP)CosS2 


V2 = ^2^2 


(7.3) 

(7.4) 


Where AP is now given by AP = OP + OA and is once again deter- 
mined from geometry. 

Then, for no gross slip, taking the cage tangential velocity due 
to roller, j, as the average of the tangential velocities 
and V 2 we get, 


Vj^g = ^{'^1 + V2) - ^^(PtQ^g) 

The cage rotation speed due to roller, j, is then given by. 


n 


eg 


,j = (^1 + V2)/Pt 


(7.5) 


(7.6) 
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ROLLER AT INNER RACEWAY 



FIGURE 9 



If is the momber of rollers loaded at both raceways, thus 

controfling the cage rotation, the expected average cage rotation 
speed is taken as the arithmetic mean of the cage speeds due to 
each of the controlling rollers as given by. 


- ( S'Scg,j'/2ieff 


(7.7) 


eg 


j=i 


7.2 ROLLER ROTATION SPEEDS 

Assuming pure rolling at the most heavily is 

inner and outer contacts, the rotational speed of a roller, :, is 
taken as the average of the two speeds imparted to this roller by 

the two raceways. 

Once again considering the contact of a roller, _j, at the 
as shown in Figure 9. The radial distance of point A, the most 
heavily loaded point, from the roller axis is given by. 


To 1 = {OA)Cos(Si-S) 


(7.8) 


For no slip at this point, the rotational speed of the roller due 
to the inner raceway contact, is then given by. 


Q 


o, 1^0,1 ~ (^I'^cg^l 


(7.9) 


Similarly, by considering the contact of the same roller at the 
outer, as shown in Figure 10, we can write. 


and 


r - = ( 0 A)Cos(S- 62 ) 
^ 0 , 2 ^ 0, 2 = (^ 2 ‘° cg ^^2 


(7.10) 

(7.11) 


Where Qq 2 rotational speed of the roller due to the 

outer ra6eway contact . 

The average of two roller speeds, due to inner > 

and outer rotation (Q^ 2 ^ equations (7 9) and T7 11) , 

gives us a good estimate of the expected roller rotational speed. 
Therefore, for the roller, j. 


°o = ^^°o,l °o,2^ 


(7.12) 


This is true for all the rollers loaded at both raceways. For 
the rollers loaded at one raceway only, the rotational speed is 
taken as the speed imparted by the contacting raceway, a 

pure rolling at the most heavily loaded point of this contact. 
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ROLLER AT OUTER RACEWAY 



INTIAL CX3NTACT ANGLE ( |3) 
OPERATING CONTACT ANGLE 

CONTACT CURVATURE CENTER 


FIGURE 10 
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8.0 RELATIVE SLIDING AT CONCENTRATED CONTACTS 


Having determined the rotational speeds of the bearing cage about 
S^ blaring axis, and of each jroller ^out its o^ax^s w^^^^ 
now estimate the magnitude and direction of expected relative 
sliding at each concentrated contact. 


8.1 INNER RACEWAY CONTACTS 

Let Oxyz be a local coordinate system with at 0 as sho^ 

in Pioure 9 The y- axis being perpendicular and into the plane 
ot ForTpoint (x,y,z), lying on the contact ellipse of 

this roller at the inner raceway, the y coordinate ^ 

as compared to x and z coordinates. This is due to the fact that 
^tact ellipse has a high ellipticity, its minor axis being 
very small compared to its maj or axis • 

Tf»t the Doint (x.z) on the contact ellipse lie on the inner 
raceway surface. The tangential velocity of this point along the 

y-axis is given by, 

= - (OP-OE) Q]_CosS 3_ + xQ]_SinS]_ (8.1) 

When the same point (x,z) is ?°^sidered lying on the 
surface, its tangential velocity along the same axis (y axis) 

would be given by. 


Vq = (OE)nQCos(Si-S) + xQQSin(Si-S) 


( 8 . 2 ) 


where OP and OE are determined from the known 

relative sliding velocity of this point (x,z) is then given by. 


^rel,l - 


(8.3) 


8.2 OUTER RACEWAY CONTACTS 

Similarly by considering a point (x,z) on the contact area of a 
roller at* the outer raceway, as shown in Figure 10, we can write 
the expressions for the tangential velocities 
lying on the outer raceway and roller surfaces as folio 

V2 = - (OP+OE) 02^03^2 xQ2SinS2 
Vq = - (OE) OqCos (IS-S 2 ) ■ xQQSin(E-S2) 

Relative sliding velocity of this point (x,z) is then given by, 

(8.3) 


■^rel,2 - ^^2 ■ ’^o 
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9.0 HEAT GENERATION IN THE BEARING 


Total heat generation in a spherical roller bearing comes from 
various friction mechanisms occurring during the bearing opera- 
tion. These include shear stressing of the EHD oil films due to 
the relative sliding motion at the concentrated contacts, viscous 
friction torque of the lubricant against the ploughing motion of 
the rolling elements through the lubricant, sliding friction at 
the cage lands and rails. In a high speed bearing of this type, 
the major contribution to the total heat generation, under normal 
operating conditions, is due to the sliding friction at the 
concentrated contacts under EHD conditions. 

In the formulation for SASHBEAN conputer program, heat generation 
due to the three mechanisms mentioned above has been considered. 
These are described in detail in the following sections. 


9.1 DUE TO RELATIVE SLIDING AT CONCENTRATED CONTACTS 

To estimate the traction force at the sliding Elastohydrodynamic 
contacts, a formulation presented by Allen et. al. [1] has been 
used. 


9.2 TRACTION COEFFICIENT UNDER EHD CONDITIONS 

As per the above model, four parameters, namely - ambient ^so- 
lute viscosity, pressure-viscosity coefficient, a lubrication 
factor (a pseudo coefficient of friction) , and a transition shear 
stress, can quantify the traction in a sliding EHD contact. 
Mathematically stated, the shear stress for a Newtonian fluid in 
a concentrated contact under EHD conditions is given by, 

T = ^^rel/^nin ^ < "^tr ^ 

= when t > t^j. and t > Mfi^ (9.2) 

Where S is the normal contact pressure at a lamina, being the 

transitional shear stress (typically 1000 psi) , and is the 

lubriccint factor - the pseudo coefficient of friction (typically 
between 0.045 and .075) . 

Having determined the shear stress (t) in the EHD film at 
each lamina contact, the total traction force for the given 
lamina contact is then evaluated by integrating the shear stress 
over the entire contact area of this lamina contact. 

With the relative sliding velocity at each lamina contact already 
known, the heat generation due to relative sliding is then given 
by the thermal equivalent of the mechanical work done against 
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this friction force. 

Rate of Work Done = Traction Force x Relative Sliding Velocity 

The thermal equivalent of this mechanical work is then taken as 
the rate of heat generation in the bearing and is given by. 

Rate of Heat Generation = Rate of Work Done/2.5933 

where the factor 2.5933 converts the mechanical work (in-lb/sec) 
to equivalent heat units (Btu/hr) . 

To estimate the sliding friction force at the concentrated con- 
tacts under lost lubrication conditions, it is assumed that the 
Newtonian viscosity relationship for the fluid shear stress no 
longer holds valid. A constant (user supplied) coefficxent of 
traction coefficient is used for evaluating the slidxng frxctxon 
forces, the mechanical work done, and the resulting heat genera- 
tion at each concentrated contact. 


9.3 SLIDING FRICTION AT CAGE GUIDING RAILS /LANDS 

The Petroff's equation [25], which provides a good approximation 
for the power loss in lightly loaded journal bearxngs, has been 
used to estimate the resisting torque of the Ixibricant present xn 
the clearances between the ring lands and the cage raxls. Fol- 
lowing assumptions have been made for this formulatxon. 


(a) Resultant cage -roller loads are very small. 

(b) Cage is properly balanced and while rotating maintains a 
uniform gap with the guide ring. 

(c) Radial gap between the ring land(s) and the cage rail(s) is 
fully flooded with the lubricating fluid. 

(d) Viscosity of the lubricating fluid in the gap does not 
change appreciably from the bulk of lubricant in bearxng 
cavities . 


The tangential friction force at the cage- ring interface is then 
given by. 

Where i=l for an inner guided cage and i=2 for an outer guided 
cage. W„_ is the total width of the cage rails in sliding contact 
with the guide ring rails, is the mean interface dxameter, 

is the angular velocity of the guiding ring, and dc^^, xs the 
diametral clearance between the guiding ring lands and the cage 
rails. The resulting heat generation due to this frxctxon power 
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(9.4) 


loss is then given by, 

= 0.1928 (D^,rf cr I Qcg"°il^ 


9.4 VISCOUS FRICTION TORQUE DUE TO LUBRICANT 

The viscous drag torque, experienced by the rollers moving 
through the lubricant flooded bearing cavities, has been estimat- 
ed using the empirical equation suggested by Harris [9] and 
Eschmann [5] for high speed bearings as presented below, 

Tfi = 1.42E-5¥{“j^N)^/^Pt^ ; > 2000 (9.5) 

Where, 

■p _ viscous friction torque of the lubricating fluid (in-lb) 

ilLi 

Y = A factor based on the bearing type and its lubrication 
system. As per ref. [9], Y = 5 for spherical roller 
bearings under oil bath lubrication. 

The heat generation in the lubricant due to its viscous friction 
torque is then given by the thermal equivalent of the mechanical 
work done against this viscous drag. This is given by, 

= 0.0404 (N.Tf]_) 0-6) 

where the factor of 0.0404 converts the rate of mechanical work 
(in- lb. /min) to equivalent rate of heat generation (Btu/hr) . 
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10.0 THERMAL ANALYSIS 


With the amount of heat generation at different points in the 
bearing now known, thermal models are then set up to perform a 
steady-state and transient heat transfer analysxs. These ana y- 
ses can predict the expected steady- state and txme- transient 
temperature maps of various points in a bearing system respec- 
tively. Some detailed information about the bearing application, 
including its supporting system, dimensional and material data, 
lubricant properties, and lubrication system, are required to 
prepare the model data for thermal analysis . 

A steady- state condition would be one when a certain l^ricant 
flow rate is maintained through the bearing system. After an 
initial "warm up" time, the temperatures at different points 
the bearing and its supporting mechanical _ system would have 
stabilized and there is no further variation in temperatures at 
any point throughout the system. 

When a "stabilized" system is subjected to a sudden change in 
environment and its steady- state condition disturbed, the system 
temperatures would become time variant. One such change o 
environment could be the loss of Ixibricant flow through the 
bearing system. A transient heat transfer analysis would be 
called for in such a situation to predict the time- temperature 
history of various points in the system, "^is information ^y 
then be used to predict bearing's time -to -failure in the event of 
lubrication failure. 


10 . 1 METHOD OF HEAT TRANSFER ANALYSIS 

A method known as "Lumped- Heat -Capacity" method has 
modeling and analyzing the system for steady- state and transient 
temperatures. In this method, the whole system is considered 
made up of small elements, with the entire thermal capacity of 
each element "Itimped" at its center and assuming a uniform tem- 
perature distribution throughout the volume of the element. 

In other words, the internal resistance to heat flow within an 
element is considered negligibly small compared to its externa 
resistance to heat flow from the surface to the surrounding 

in general, snaller the size of elements, the more realis 
tic these assumptions are for lumped-heat- capacity analysis. 

For ease of data preparation, the maximum number of elements for 
a single bearing system have been limited to 20 in this program. 
The model prepared for the steady- state and transient therma 
analyses of McGILL SB- 1231 and SB- 1231-1 bearinp is shown in 
Figure 11. The element and node numbers are marked on the s^e 
figure. The lubricating oil circuit is shown by the dotted lines. 
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THERMAL ANALYSIS MODEL FOR THE 
TEST BEARING AND ITS SUPPORTING SYSTEM 






10.2 fiTRADY- STATE ANALYSIS 

SASHBEAN Con^uter Program is capable of analyzing and predicting 
the steady- state temperature map of an axisymmetric mechanical 
system of any cross-section. The mechanical system is first 
approximated by an equivalent system comprised of a number of 
ellJiLts of simple geometries. Each element is then represented 
by a node point and its node munber. The tenperature at this 
node point (element) may or may not be known. Heat sources (with 
known heat generation rates) are assigned to any nodes represent 
ing the elements of heat generation within the The 

environment surrounding the system, also having heat transfe 
with the system, is also represented by one or more nodes. As 
mentioned earlier, a maximum of 20 nodes are allowed by the 
program to describe the equivalent thermal model- 

with the elements and their node points properly selected, the 
heat balance equations considering different modes of heat trans- 
fer are then set up for each node. Heat transfers by conduction, 
free and forced convection, and mass transfer are considered. 
Radiation heat transfer, being very small as compared to other 
modes in such applications, has been neglected. 

Now for a steady state condition to exist, the net flow of heat 
to a node i from its surrounding nodes j plus the heat generated 
at node 1 must be equal to zero. This heat balance at node i 
results in an algebraic equation with some nodal teinperatures as 
unknowns. As this condition would be true for all the nodes in 
the system, heat balance condition at each node results in an 
algebraic equation for that node. 

Assuming a linear relationship for free ^d forced convection 
modes, the resulting mathematical model is a system of linear 
algebraic equations with the unknown nodal temperatures the 
unknown variables. The analysis then reduces to soO^ving this 
system of algebraic equations for the un]oiown variables. 

Consider the heat energy flowing into node 

nodes j . If N is the total number of nodes in the model , then 
for the steady state condition to exist. 


Net heat flowing into node i = o 


or. 


Qi 


n 
+ E 

j 



+ f4(Tj-Ti) 


0 ( 10 . 1 ) 


Where f^, 
transfer 
elements . 
provides 


f- f^ and tn are the coefficients in respective heat 
ecfuations for various modes of heat transfer between the 
Volume I of this documentation, the USER'S GUIDE, 
more details on various modes and their applicable 
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( 10 . 2 ) 


equations . Rewriting equation (10.1) we get, 

Qj, = E Fj^j(Tj-Tj_) = 0 ; i = 1, n 

where , 

^ij = ^1 ^2 ^3 ^4 

Equation (10.2) represents a system of linear algebraic equations 
for 1 s i s N. Depending upon the actual ntimber of surrounding 
nodes j interacting with node i, each equation may or may not 
have all the unknown variables. 

For a node (say node j) with a given (known) temperature, the 
heat balance equation for this node is replaced by the equality 
ip . _ T • and the unknown variable Tj sxabstituted for the 
kAown,^^*^^^^^^, in the remaining equations of the system. 

The Gauss -Jordan numerical scheme, with partial pivoting, has 
been deployed to invert the coefficient matrix for the solution 
of this system. Reference [3] provides a detailed discussion on 
this and other numerical methods. 


10.3 TRANSIENT ANALYSIS 

The formulation for computing the steady- state temperatures was 
based on the condition that the net ener^ transfer into any 
node, i, from its surrounding nodes, j, is zero. Therefore, no 
further change of element temperatures then took place at this 
condition. 

For the transient formulation, on the other hand, a net energy 
transfer to the node i from its surrounding nodes j takes place 
resulting in an increase in the internal energy of the element i . 
Precluding any phase changes, this increase in the internal 
energy results in a temperature rise for the element i. As each 
element voliime is considered to have "lumped" capacity at its 
node point, the interaction of all the elements thus determines 
the behavior of the complete system. 

For the mathematical formulation of the transient problem, the 
net heat energy into a node i from its surrounding node(s) j in a 
very small interval of time (dt) is equated to the energy re- 
quired to raise the temperature of the element i to a new value. 
Heat transfer by conduction, free and forced convections are 
considered. To simulate the lost-lxibricant condition no heat 
transport by mass transport of the lubricant is now available. 

The lubricant nodes in the steady state model are replaced by air 
nodes having a forced convective heat transfer with the bearing 
elements . 

The resulting equation for the node, i, is a linear differential 
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eauation of first order with a known initial value. The ele- 
ment's steady-state temperature, estimated prxor to thxs analy- 
siS/ is used as the initial value at tiine t = o. 

If the internal energy of an element, i, is e^ressed in te^ of 
its specific heat and temperature, its rate of change with time^ 
t is equal to the net heat energy gaxned by the node, x, xn the 
sxasill interval of time, dt. Therefore we can wrxte that. 


dEi dT^ 

q (CiPiVi) — 

dt dt 


(10.4) 


Where c is the specific heat, p the mass densxty, and v the 
volume of the element, i, material. Usxng an . 

difference approximation we can then transform _ the dxfferentxal 
equation (10.4) into a finite difference equatxon as follows. 


> ^ p+1 _iji . p 


Delta (Ej^) T^^ 

q = — = (Cdp^v^) 

Delta (t) Delta (t) 


(10.5) 


Where tP and tP'^^ are the temperatures of element, i, at time t 
and t+Delta(t) respectively. 

The total heat gain, q, of this element, bv%!rioSs'"'^^''^ 

elements, j, in the same txme xnterval. Delta (t), by varxous 

modes of heat transfer, is also given by. 


q = Qi 1 E Fij (TjP-TiP) 


( 10 . 6 ) 


Therefore by equating (10.5) and (10.6) and solving for Tj^P we 
get. 

Delta (t) Delta (t) 

T-P"*"^ = (Q4 + T FjjT^P)^ + [1 - • S 

1 ^ j ^ ^ Ci Ci D 

Where, = c^p^Vj^ is the thermal capacity of the element i. 

If T-P is known at time t, then Tj^P”*"^ at time t+Delta(t) can be 
detemined from equation (10.7) for the element x. Thxs process 
iSrepeated for each node in the model and thxs txme-marchxng 
solution continues till the desired t^erature or 
reached. As mentioned before, the nodal steady- state temper 
tures are used as the initial values at time t = o. 
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APPENDIX - A 

VOLUME OF A FULLY CROWNED SPHERICAL ROLLER 


VoliJine of a spherical roller, required in the estimation of 
centrifugal force and gyroscopic moment acting on a roller, is 
determined using simple integral calculus as o ows. 

Let OXYZ be a coordinate frame of reference with its origin at 
the roller center as shown in Figure 12. The coordinates of the_ 
roller crown radius center (O') are then given by Xq - 0 and y^ 

- (R -D/2). Therefore, the profileof the roller crown can be 
represented by the following equation of a circle. 


or. 


(x-Xq) 2 + (y-yo)^ = 

y = (Rq^-x^)^ + : 


(A.l) 
(A. 2) 


Now by considering the volume of a roller slice of differential 
Sictaess, dx. at a distance x from the roller 

Shown in Figure 12) and integrating it for the full roller lengtn 
we get. 


V = 


dV = IT 


y2 dX 


(A. 3) 


By substituting for y = y(x) from equation 

(A. 3) and integrating w.r.t x we get an expression for the volume 
of a fully crowned spherical roller. 
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APPENDIX - B 

MOMENTS OF INERTIA OF A SPHERICAL ROLLER 


Mass Moments Of Inertia (MMOI) of the roller both its 

longitudinal and transverse axes, are re^xred for ^e 
tions of gyroscopic moment acting on each roller. The e^r 
sions for these have also been derived by ^ 

slice of differential thickness and integrating over the full 
length of a roller, as summarized below. 

By referring to Figure 12 and considering the roller slice as 
shown, the moments of inertia of the slice about the X, Y, and 
coces are given by, 

(B.l) 
(B.2) 
(B.3) 


dlv = i^y^dm 


dly = dly. + x^dm = (^2^x2) dm 

din; = dl 2 . + x^dm = (",^2+x2)dm 


Where divi and dl^, are the moments of inertia of the disk about 
a local coordinati axes X'Y'Z’, with origin at the disk center 
and parallel to the global axes OXYZ as shown in Figure 12. dm 
is the mass of the differential disk and is given by, 

dm = pdv = pTry^dx (B.4) 

p being is the mass density of the roller^mategial . , 

spherical profile is described by y = (Rq -Xg) ^o 
in equation (A. 2) of Appendix-A with yQ = -(Ro-D/2). 

By substituting equation (B.4) into equations 

(B.3) and performing the integration w.r.t. x for x ry ^ ^ 

-L/2 to +L/2, we get the expressions for roller's mass moments of 
inertia about its longitudinal axis (I^) and transverse axes (ly, 
I ) . Due to its symmetry about the X-axis, we get ly ij;- 
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APPENDIX - C 


INITIAL ESTIMATION OP INTERNAL ROTATIONAL SPEEDS 

To estimate the magnitudes of the roller dyiiamic loads, namely 
the centrifugal force and gyroscopic moment, cin initial estimate 
of the cage rotational speed cdDout the bearing axis and each 
roller's rotational speed about its own longitudinal axis was 
first made under the following assuit5>tions ; 

(a) The operating contact angle of each roller is the initial 
design contact angle of the bearing. In other words, the contact 
ellipse at each concentrated contact is centered about the roller 
center. 

(b) No gross slip occurs at any concentrated contact. 

The tangential velocity of the cage is then given by the mean of 
the tangential velocities of contact center points at the inner 
and outer raceways. Thus, 

Q^g = ^[Qid-e) + Q2(1+0)] (C.l) 

or, 

N^g = MENj^d-e) + N2d+e)] (c.2) 

The rotational speed of each roller about its own axis, required 
for estimating the gyroscopic moment on each roller, is then 
given by, 

Nq = (Pt/2D) (1-0) (1+0) (N 2 -N 1 ) (C.2) 
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appendix - D 

ROLLER GYROSCOPIC MOTION ANALYSIS 


Tn anoular contact spherical roller bearings, the motion of each 
^ol“f ifve?? to that of a gyroscope ^ 

rotor, the rSler is rotating (gyroscopic spin) S^bit- 

aeometric axis, unparrallel to the bearing ^is, also o^bit 

ing (gyroiscopic precession) about the bearing axis along with 
the caae The gyroscopic moments seen by these rollers, 
neolec?Ia speed and/or low contact angle analyses become 

^?te sfgnScant in high speed applications and consrderably 
affects the roller deflections and load distributions. 

The formulation for the SASHBEAN Computer Program takes into full 
consideration this moment loading of the rollers. The a p 
tions S ?or the estimation of the magnitude of this moment 
acting on each roller are listed below. 

(a) The cage has a constant angular velocity of rotation about 
its (or bearing) axis as determined in Appendix C. 

(b) Each roller has the same constant angular velocity about its 
longitudinal axis as determined in Appendix C. 

(c) The roller misalignment (pitching) angle at 

location is negligibly small when compared to the design 
contact angle of the bearing. 

(d) The operating contact angle of each roller is the same as 
the design contact angle of the bearing. 

Consider a roller at any angular location (0) 

Fioure 13 Let OXYZ be a fixed frame of reference with its 
orS?n, o; at the intersection of bearing and roller axes. Let 
o’xyz be a rotating system of axes, attached to the ^ _ 

origin at the roller mass center as shown in the 

considering the roller as a gyroscopic rotor we can see that fo 
this gyroscope, 

Nutation Angle = Roller -Raceway Contact Angle = S 

Precession Angular Velocity = Roller Orbital Velocity = Q^rg 

Spin Angular Velocity = Roller Rotational Velocity = Qq 

Bv considering the special case of the gyroscopic motion where S, 
^ , and remain constant, the couple (moment) retired on the 
roller (the rotor) to sustain this motion is given by. 


GM = [Iz ^ 


(D.l) 
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This couple, applied to the roller by the raceways, is about an 
axis perpendicular to the precession and spin axes of the ^ro- 
scope. The moment vector is thus directed along the positive y- 
axis (perpendicular and pointing into the plane of the paper) . 

The reaction couple of the roller, resisted by the raceways, is 
ec[ual in magnitude and opposite in direction. For a detailed 
discussion on the topic and derivation of the above equation, the 
reader is referred to reference [2] . 



FIGURE 13 
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APPENDIX - E 


LOAD -DEFORMATION RELATIONSHIP AND MATERIAL FACTORS 

For a cylindrical body of finite length, when pressed onto a 
plane sLfaced body of infinite extension, the no^l approach 
between the axis of the cylinder and a distant poxnt ^he 
supporting body, first presented by Palmgren et. al. [22], 
approximately given by , 


or 


8 = 4.36E-7[Je2*'^Q-^/L*®] 
5 = 4.36E-7[iE2*'^q*^L*^] 


(E.l) 

(E.2) 


where iE is a factor based on the materials of the two contacting 
bodies and Q = qL is the total normal force pressing the t 
bodies together. The material factor JE is given by. 


E 




(E,+E2)-i^'^^ 
643E+7-^^ 


(E.3) 


E1E2 


where, E^^ = E^^/d-e^^) and E 2 = E 2 /(l-e 2 ^)- 


(a) RELATIONSHIP FOR STEEL ON STEEL CONTACT : 

Let = 29.0E6 psi, E 2 = 29.0E6 psi, e^^ = 0.30, and e 2 
we get = 1.0. Equation (E.2) then gives us, 

5 = 4.36E-7cr^.L*^ 


= 0.30 


(E.4) 


(b) relationship for ceramic on steel contact : 

Let E^ = 43.0E6 psi, E 2 = 29.0E6 psi, e^ = 0.27, and e 2 = 0.30 
We get = 0.872. Equation (E.2) then gives us, 

6 = 3.80E-7.q*^.L*^ 
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